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Abstract 

Heteroepitaxial self-assembled quantum dots (SAQDs) will allow breakthroughs in elec- 
tronics and optoelectronics. SAQDs are a result of Stranski-Krastanow growth whereby a 
growing planar film becomes unstable after an initial wetting layer is formed. Common sys- 
tems are Gea;Sii_a;/Si and Ina;Gai_2;As/GaAs. For applications, SAQD arrays need to be or- 
dered. The role of crystal anisotropy, random initial conditions and thermal fluctuations in 
influencing SAQD order during early stages of SAQD formation is studied through a sim- 
ple stochastic model of surface diffusion. Surface diffusion is analyzed through a linear and 
perturbatively nonhnear analysis. The role of crystal anisotropy in enhancing SAQD order is 
elucidated. It is also found that SAQD order is enhanced when the deposited film is allowed 
to evolve at heights near the critical wetting surface height that marks the onset of non-planar 
film growth. 

1 Introduction 

Heteroepitaxial self-assembled quantum dots (SAQDs) represent an important step in the advance- 
ment of semiconductor fabrication at the nanoscale that will allow breakthroughs in optoelectron- 
ics and electronics. [1-12] SAQDs are the result of a transition from 2D growth to 3D growth 
in strained epitaxial films such as Sia:Gei_a./Si and In2.Gai_xAs/GaAs. This process is known as 
Stranski-Krastanow growth or Volmer- Webber growth. [1,13-15]. In applications, orderof SAQDs 
is a key factor. There are two types of order, spatial and size. Spatial order refers to the regularity 
of SAQD dot placement, and it is necessary for nano-circuitry applications. Size order refers to 
the uniformity of SAQD size which determines the voltage and/or energy level quantization of 
SAQDs. It is reasonable to expect that these type of order are linked, and it is important to under- 
stand the factors that determine SAQD order. Further understanding should help in the design and 
simulation of both spontaneous "bottom up" self-assembly and directed or guided self-assembly to 
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enhance SAQD order. [16-23] Order prediction inherently involves the modeling of stochastic pro- 
cesses. Recently, SAQD order has been modeled using a deterministic model with stochastic initial 
conditions in the linear approximation. [24, 25] This model was based on previous nucleationless 
models of SAQD formation. [26-29] Here, this method of modeling SAQD order is improved by 
incorporating stochastic thermal fluctuations in the surface diffusion. Thus, the previously deter- 
ministic governing equations become stochastic. The final order predictions are qualitatively the 
same as for the previous linear model, but they are quantitatively different. Additionally, prelimi- 
nary non-linear modeling results are presented. One non-linear model approximates a ID surface, 
but incorporates the stochastic thermal fluctuations. The second non-linear model is of a 2D sur- 
face, but it is only implemented as a deterministic model at present. 

In the previous work using a linear deterministic model with stochastic initial conditions [24, 
25], it was found that peaks in the linear dispersion relation can be used to predict and explain 
order. It was also found that only anisotropic models give rise to dispersion relations with dis- 
crete peaks, thus explaining why elastic anisotropy contributes to SAQD order as previously re- 
ported. [28, 30-32] The dispersion relation was then used to generate a spectrum function in the 
hnear approximation, and the spectrum function in turn could be used to define and predict two 
correlation lengths that grow as the square root of time. These correlation lengths were identified 
as the key quantities describing SAQD order. Using equations for these correlation lengths, it was 
found that growth of SAQDs with an average film height near the critical 2D-3D transition height 
might enhance order, although practical limitation of producing ordered arrays of SAQDs were 
also noted. Although the incorporation of a wetting potential is possibly controversial, it appears 
to produce the correct phenomenology, and it may possibly be a mathematical surrogate for more 
complicated processes such as stabilization by intermixing. [33] See refs. [25-29, 34] for further 
discussion. In ref. [25], it was also shown that the modeling/order prediction method could easily 
be applied to a large class of models, but the simplest model that produced Stranski-Krastanow 
growth was used as an example. Additionally, various mathematical issues such as convergence 
and precise definitions of the correlation functions as either spatial averages or ensemble averages 
was treated. Readers interested in these more technical details are referred to ref. [25] . 

The new result presented here is mainly the mentioned incorporation of thermal fluctuation 
to seed quantum dot formation, as opposed to the somewhat artificial assumption of a random 
roughness initial condition that is chosen more or less arbitrarily. One product of the present work 
is a formula to choose this initial roughness to give a nearly equivalent disordering effect as thermal 
fluctuations; however, a deterministic model will never be a true substitute for a stochastic one. 
The outcomes of the stochastic model are qualitatively similar to the previous deterministic model, 
but quantitatively distinct. In addition to the stochastic linear model of SAQD order, preliminary 
results of non-linear models are presented. These models appear to corroborate the linear model 
predictions but also give a more complete picture of the time evolution of SAQD order. The 
current model predicts that order will be fairly poor under most growth conditions. This seems to 
be in agreement with most experiments, for example refs. [14,35-38]. The basic phenomenology 
appears to be more or less in agreement with observations; however, more quantitative reporting 
of experimentally observed order would facilitate future comparisons. 

The rest of this article is organized as follows. Section |2] presents the stochastic governing 
equations and physical causes of SAQD formation. Section [3] presents the linearization of the 
model presented in Sec. |2] along with the extraction of order predictions and application to growth 
near the critical film height using parameters appropriate to Ge/Si SAQDs. Section |4] presents 
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preliminary non-linear modeling results. Section [5] presents the conclusion. Finally, Appendix [A] 
presents the derivation of the time evolution equation of the spectrum function. 



2 Physical Model 

The formation and growth of SAQDs is modeled in a fashion similar to refs. [18,26-28,39]. The 
film surface is described by the film height as a function of the lateral position, 7Y(x). The film 
height evolves via surface diffusion that is driven by a diffusion potential, /x(x). The film surface 
grows with a velocity normal to its surface that is given by 

f „ X = nJx' 



dt 

VVsfi{^) + ^/2QVkbTr]{^,t) 
••• + n^(x)Q, (1) 

where V is the surface diffusivity; nz{x) is the z— component of the surface normal vector, n(x); 
Vs is the surface gradient; 'Vs- is the surface divergence; Q is the flux of new material onto 
the surface; and y/2Q,VkhTrj{x, t) is the fluctuation of the surface diffusion. Note that the sur- 
face diffusivity is assumed to be a scalar; thus, it is isotropic. A limited discussion of diffu- 
sional anisotropy appears in ref. [25], and full development is in progress. The surface diffu- 
sion fluctuation is chosen to give a steady state that is consistent with the Gibbs distribution for 
a quadratic potential. [40] In ref. [40], there is a slope-dependent intensity factor, but here that 
factor is neglected for simplicity and because it has no effect to linear order. r]{x,t) is a time 
fluctuating white noise (or the derivative of a Brownian process) [41,42] so that it has zero mean 
(r7(x, t)) = 0, where (...) denotes the ensemble average, and it has a sharply peaked correlation 
function, (r7(x, t)r}('x', t'))_= I5^(x — x')5{t — t'). I is the rank 2 identity matrix, and 5{x) is the 
Dirac Delta function. Eq.[l]is interpreted as an Ito stochastic differential equation}^ 

The diffusion potential is derived from the total free energy. The details of the derivation are 
covered in ref. [25], and only the most important points are reviewed here. The total free energy 
is assumed to have two parts, elastic energy and a term that is a combined surface energy and 
substrate wetting energy JF = jFgiast. + J^sw The second part Tsw is an integral over the horizontal 
coordinate x of an areal energy density. 



rs^= I c^2xF,^(7^(x),v7^(x)). 

J X— plane 



The areal energy density, Fsw, is in turn a function of the film height, 7Y(x) and the film height 
gradient, V7-^(x). From this total free energy, one can find the diffusion potential ;u(x) by taking 
the variational derivative with respect to film height and multiplying by the atomic volume, Vt, 

Mx) = ^]Kx) + i^^°)(x)-V■F(°J)(x)]. (2) 



'Because the noise term is additive and not multiplicative, it does not matter whether Eq.[l]is interpreted as an Ito 
stochastic equation or Stratonovich stochastic equation. [41,42] 
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co'(x) is the elastic energy density at the film surface. Fsw^'' indicates the derivative with 
respect to H and the n"^ derivative with respect to V7Y. Fsw^\x.) = d'H(yi)Fsw ('^(x), V?^(x)) and 



each vector component of Fi'i}''(x) is Fs^^''(x) = i9[vw(x)]^-^su> ('^(x), V7i(x)). This diffusion 
potential (Eq. [2]) is a general form for any surface diffusion model that incorporates the non-local 
elastic energy density and a local areal energy density such as a surface energy (even one with 
orientation dependence/faceting [27]) and a wetting energy. [27-29,34]]^ In refs. [24,25] a simple 
model is analyzed that includes elastic anisotropy, a constant surface energy density, 7, and a 

substrate wetting energy density, W{7i). For this simple model, Fsw = [l + (^"^(x))^] 7 + 
W (7i(x)), and the resulting diffusion potential is 

/i(x) = n [^(x) - 7«:(x) + W'iUi-x))] , (3) 

where k(x) is the total curvature, and W'{T-C) is just the derivative of the wetting potential. A more 
extensive discussion of different possibilities for Fsw is discussed in ref. [25]. 

3 Linear Stochastic Model 

Stochastic terms that fluctuate in time lead to stochastic differential equation that are often difficult 
so solve with either analytic techniques or numerical simulation. [41,42] Linear stochastic differ- 
ential equations, however, are much easier to solve. In fact, their solution is not very different 
from the solution of linear deterministic (or ordinary) differential equations with stochastic initial 
conditions. The linear model is naturally more approximate than the non-linear model, but it rep- 
resents an important first step, and its solution can facilitate the development and interpretation of 
non-linear models. 

To model the development of SAQD order, the growth dynamics are linearized producing a 



^(01)/ 



linear dispersion relation (Sec. 3.1 ). Then, the spectrum function is calculated based on the gov- 



erning linear equations and the dispersion relation (Sec. |3]2). The expression for the spectrum 



function is then applied to the simple diffusion potential (Eq. [3]) for a (100) surface of a cubic 



crystal (Sec. 3.3 1. Application of this method generally (to other surfaces or crystals) is outlined in 
ref. [25]. In this part of the calculation, crystal anisotropy can play an important role in the diffu- 
sion dynamics and development of SAQD order. [24,25,28,43] For simplicity, it is assumed that 
only elasticity has a strong anisotropic effect. A more detailed analysis of other anisotropic effects 
can be very cumbersome. [25] Using the specific dispersion relation, formulas for the correlation 
lengths that quantify SAQD order and the real-space correlation function or derived. Finally, the 
correlation function and the correlation lengths are applied to a numerical example of Ge dots on 
a Si substrate. In this example, order predictions and dependence of order on average film height 
is compared with previous deterministic models. 

3.1 Linearized IModel 

Eqs.[T]and[2]are linearized about the average film height (denoted H) for the case of zero deposition 
rate (Q = 0). Thus, the following analysis would correspond to a fast deposition and then an 



It is possible that the wetting potential is simply an approximation to the stabilizing effect of intermixing. [33] 
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Figure 1: Evolving film surface. Total height is average (7i) plus fluctuations h. 

anneal. Other growth cases such as constant deposition rate can be analyzed in a similar fashion, 
but they are beyond the scope of the present work. Following refs. [29,44], the total film height is 
the average film height plus small fluctuations (Fig.[T]), 

H(x,t) =7i + /i(x,t). 

Due to translational invariance of the governing equations, the Fourier components of /i(x, t) 
evolve independently in the linear model. Also, the non-local nature of the elastic energy makes 
calculations using Fourier components (spectral methods) easier than using /i(x, t). Fourier trans- 
forms use the convention, /(x) = J (i^ke*'^'^/k and /k = (27r)~^ J (i^xe~*'^'^/(x). /ik is the 
Fourier transform of /i(x), where k is the corresponding wave vector. 

The linearized diffusion potential is calculated following ref. [25]. Linearizing the surface- 
wetting part of the diffusion potential, Eq.|2]and taking the Fourier transform, one gets [24] 

/is.,ii„,k = n + k ■ F(°J) ■ k) 

where the Fsw^^ terms are the derivatives of Fswi^H, V7Y), evaluated for a perfectly flat surface of 
height H. They are constants in the following analysis because they depend only on the average 
film height H. The first superscript indicates the derivative of with respect to H. The 
second index indicates the n"^ derivative with respect to 'VH. Evaluated for a perfectly flat surface 
of height n, Fir"^ = d^di-y^Fs^ {n, vn)\ -)i^=-fi^ vw=o • The elastic energy density at the film sur- 
face is calculated as in refs. [25,45,46], where the bimaterial (film -i- substrate) is approximated as 
an elastically homogeneous material to simplify calculations]^ The resulting elastic energy density 
to linear order is cona^u = —Sg^khi^ so that the elastic energy is proportional to the wavenumber 
k = ||k|| and the Fourier component /i^, and it has a prefactor that depends on the wave vector 
direction, Oy^. Thus, the total linearized diffusion potential in reciprocal space is 

/iiin,k = n (^-Se,k + + k ■ FfJ) ■ k) /ik. (4) 

Linearizing the dynamic evolution, Eq. [T| and plugging in Q = and /iiin,kj 

dM) = (Tk/ik(t) + ^2nVhT [zk • 77k(t)] ; (5) 
c^k = -fc^Pl)(-^,,A: + i^^°) + k.F(°J).k), (6) 

■'Ref. [44] treats a bimaterial, but an elastically isotropic one. 
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where rj^it) is the Fourier transform of r7(x, t). It has zero ensemble mean, {riy_(t)) = 0), and a 
sharply peaked two-point correlation function, (T7k(^)^k'(^')*) = (27r)"^5^(k - k')5{t - t')|^The 
growth rate of each Fourier component, a^, is dubbed the dispersion relation. 



3.2 Spectrum Function 

Eqs. |5]and[6]can be solved as a system of uncoupled linear stochastic ordinary differential equa- 
tions with constant coefficients [41,42] because the Fourier components, /ik, evolve independently 
to linear order. One could assume that there are both stochastic initial conditions and thermal 
fluctuations; however, the purpose here is to analyze the impact of just the thermal fluctuations on 
order. It is assumed the film is perfectly flat at t = 0, and that the instability is seeded by just the 
thermal noise. Thus, initially, /ik(0) = for all k, and to linear order, the ensemble average film 
height fluctuation remains zero for all time. However, the spectrum function, Ck(t) provides the 
lowest order non-trivial statistical description of film height fluctuations, and it is used to predict 
the order of SAQD arrays in a fashion similar to refs. [24, 25]. By taking the inverse Fourier trans- 



form of the spectrum function, one can predict the real-space correlation function (Sec. 3.3.1 1. A 
more complete picture of the interrelations between the spectrum function, the real- space correla- 
tion functions and other correlation functions is presented in ref . [25] . 
Taking the ensemble average of Eq. [5} 

dt {K{t)) = ak {K{t)) + ^mmr [zk • (r7k(t))] . 



The surface diffusion thermal fluctuation is mean-zero (Sec. 3.1 1, and the initial surface height 
fluctuation is mean-zero, (/ik(0)) = 0; thus, {h^it)) = for all time. 

Starting from the linearized governing equation and initial conditions, an evolution equation 
for the spectrum function can be derived that is both linear and deterministic (Appendix [A]), 

dtC^it) = 2akCk(t) + {2nVhT) . (7) 

Using the initial condition that Ck(0) = 0, 

The spectrum function Ck(t) is the average value one would expect if one extracts from a simula- 
tion or experiment the film height power spectrum, (27r)^ \ h]i{t)\^ /Area ^ Ck(t). [25] 



3.3 Application to (100) surfaces 

The spectrum function time dependence, Eq.[8} is now applied to a (100) surfaces of cubic crystals 
using the simple diffusion potential, Eq. |3| Anisotropy plays an important role in order develop- 
ment, and for simplicity only elastic anisotropy is included. From this analysis, the two correlation 
lengths are found, and then the correlation function. Finally, a numerical example of Ge dots on 
a Si substrate is presented. The dependence of order on film height is investigated and compared 
and contrasted with the similar dependence from the previous deterministic model. [25] 

"^Note that because ri{x,t) is real, J7k(t) = ri_-^{t)*, where "*" indicates complex conjugate. 
Thus,(r7k(t)r7k'(t')> = i2Tr)-^6^{li + k')<5(t - t'), as well. 
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3.3.1 Spectrum and Correlation function 

If (Tk is peaked at wave vectors, k„, corresponding to some reciprocal lattice vectors, then a 
quasiperiodic arrangement of SAQDs can form during the initial stages of growth. [24, 25] This 
quasiperiodicity is demonstrated by applying the linearized simple diffusion potential, Eq.|3| along 
with elastic anisotropy cj(x) to Ge deposited on Si with a (100) substrate surface. For a (100) 
surface of a crystal with cubic symmetry, dJun ^ = —£o° (l — e^i sin^(26'k)) kh^ is a very good 
fit to a full elasticity calculation, where Eqo is the elastic energy prefactor for 9^^ = 0°, and 
eA = {£o° — £45°) /£o° is an elastic anisotropy factor. [25] The resulting linear diffusion poten- 
tial in reciprocal space is [25] 

/xiin,k = ^ [Soo (1 - eA sm\2e^)) k + jk' + W"{n)] h^, 

where 7 is the surface energy density, and W"{H) is the second derivative of the wetting potential. 
One can see that this is a special case of Eq. |4} 
The corresponding dispersion relation is 

(Tk = Vnk'^ [So. (1 - eA sin2(2^k)) k - - W'iH)] , 

assuming that diffusivity is isotropic as in refs. [24,25]. 

From this dispersion relation, characteristic lengths and times can be found along with details 
of the early film evolution behavior. A characteristic wavenumber and time can be defined, k^ = 
Sqo/'J and = 'j^ (V^ISqo) . Also, the strength of the wetting term W"{H) can be expressed as a 
dimensionless variable, P = 'jW" (H) / £qo . A detailed analysis [24,25,27,29] shows that a large 
value of W"{H) such that /3 > 1/4 stabilizes a flat film to linear order in /ik, while a small value 
of W"(H) such that P < 1/4 is insufficient to stabilize a flat film for all possible fluctuations, h^,. 
For sufficiently small /3, cTk has 4 peaks along the four (100) directions, corresponding to 6'k = 
0°, 90°, 180° and 270° and k = aok^. = (3 + V9 - 32/?) /8 is a convenient dimensionless 
quantity. Thus, the four peaks occur at ki = aokci, k2 = aokd, ks = — ao^ci and k4 = — ao/cj. 
Expanding in the vicinity of peak n, 

ak^ao- ^(ri\{ki\ - aokcf - ^(r±kl, (9) 

where /cy is the component of k parallel to k„, and k± is the component of k perpendicular to k„. 
ao = \t~^al (tto — 2/3), a\i = t~^k~'^ (3ao — 4/3), and a± = 8eAaot~^k~'^. 

Eq. |9]is used to find an approximate expression for the spectrum function Ck(t). o"k appears 
inside an exponential; thus, for sufficiently large values aot, the exponential term in the vicinity 
of the peaks will dominate over all other contributions to the spectrum function. Thus, Ck(t) will 
have the approximate form of four gaussians each centered around the four peak locations, k„. For 
sufficiently narrow gaussians, the prefactor can be approximated by its value at the peak. Thus, 

(2vr)Vo 



^e"5^n^"~''"^ -^^x''^, (10) 



n=l 
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Figure 2: (a) Film height and (b) spectrum function of the = l-lHc simulation discussed in 



Sec. 4.2 at the end time, t = 255tc- The drawn circle in (a) has a radius equal to L± calculated 



from the spectrum function (b). 



where 



L|| = ^2^ = fc;V(6«o-8/9)(t/g, 
L± = y/2a±t = fc-Vl6eao(t/tc). 



(11) 
(12) 



L|| and L± are the two correlation lengths that arose from models with deterministic evolution 
and stochastic initial conditions. They are measures of how spatially ordered an array of SAQDs 
is. The distance over which one can expect an array of SAQDs to appear periodic is about twice 
the smaller of the two correlation lengths, usually L±. [24,25] Fig. |2] shows an example of a film 
surface with the correlation length indicated. The approximate spectrum function, Eq. [10} is only 
valid when aokcL\\ 1, and aokcL± ^ 1. Of course, when this is not the case, order will be 



very poor. Thus, Eqs. 10-12 are useful for quantifying order when it is good, and they are able to 



indicate when order is poor. 

The spectrum function, Eq. 10 is very similar to the spectrum function for the deterministic 
case with stochastic initial conditions characterized by a noise amplitude A^. [24,25] If the noise 
amplitude is set to be 

= vnhT {aoKf /ao, 

then the two cases are equivalent to linear order, when one performs these similar expansions. Of- 
ten, one uses deterministic evolution equations with stochastic initial conditions as approximations 
to stochastic evolution equations. By performing a suitable linear analysis as done here, perhaps 
one can find an appropriate initial condition for such approximations. Note that A^ has dimensions 
of [length]"^, and the size of fluctuations in a discretization procedure changes with the discretiza- 
tion length scale. The spectral methods used here handle this problem fairly easily as one can 
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coarse-grain a model by simply discarding fast oscillating noise components. A spatial discretiza- 
tion such as finite differencing or the finite element method makes quantitative implementation of 
white noise more complicated. 

As with the deterministic model [24, 25], one can take the inverse Fourier transform of the 
spectrum function to obtain the real-space correlation function, 

C7(Ax) = ^^M;^e-o*... (13) 
7rcroL||L_L 



■ ■ ■ X 



e 2V 'II " ' ^' cos[aQkci\x 

•■■ + e 2V /X cos[aokci\y) 

The correlation function, C(Ax), is a good predictor of the autocorrelation when the sampled or 
simulated area is very large. [25] 

3.3.2 Numerical Example and film height dependence 

In ref. [25], it was found that for reasonably soft wetting potentials, there can be some enhancement 
to spatial order when annealing takes place for films with heights, 7Y, that are only just above 
critical film height for unstable 3D growth. This finding was based on an assumption that the order 
that develops during the initial stages of growth is a meaningful order estimate. This assumption 
is justified to an extent by published numerical simulations [28,29,47-49] and is further justified 
by initial non-linear modeling results in Sec. |4j In ref. [25], the correlation lengths were found 
using parameters appropriate to Ge deposited on Si. A condition for the end of the linear evolution 
regime was taken to be when the r.m.s. film height fluctuation exceeded the atomic scale, the height 
of one monolayer. The r.m.s. height fluctuation is just K,ra.&. = \C (0)]^^^, using Eq.[l3 The time 
at which this condition was satisfied, tiarge^ was recorded, and the smaller correlation length, Lx_ 
was calculated for this time. These predicted tiarge values and the number of correlated dots in a 
row were graphed vs. the dimensionless wetting potential strength for /3 = . . . 0.25. It was found 
that the calculated time tiarge and the calculated correlation length diverge as (3 0.25. 

The same procedure is now followed for the present model for Ge on Si with temperature 
T = 600° C. All values are the same as for the calculations in ref. [25]. The results are graphed in 
Fig.|3j When compared with the results from the deterministic model with stochastic initial condi- 
tions [25], one finds that the effect of thermal fluctuations in the surface diffusion are actually more 
disruptive to order than assuming an initial surface with atomic scale roughness. The qualitative 
trends are the same,however, and the divergence in correlation length as /3 ^ 0.25 is observed. As 
discussed in ref. [25], one should take care interpreting this result, and there are of course practical 
hmitations. The order enhancing effect of near critical growth has not been experimentally ob- 
served (or looked for), and there may be practical limitation to implementing near critical growth 
as a method to enhance order such as the requirement for precise deposition control. 



4 Perturbatively Non-Linear Models 



The order estimates presented in refs. [24, 25] and Sec. 3.3.2 are based on the order that develops 



before fluctuations become large. The significance of these calculations is based on the following 
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and vs. f, (Ge/Si) 




Figure 3: tiaige/tc and L±/Lo vs. the dimensionless wetting parameter P for Ge/Si at T = 600°C 
as discussed in Sec. 13.3.21 



observations: 



1. Order increases during the linear stage of growth as t^/^ (Eqs. 
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and 



12). 



2. Order does not increase forever. If it did, growing perfectly ordered arrays of dots would be 
trivial. Also, qualitative analysis of numerical simulations bears this out. [28,29,47-49] 



It is, of course, worthwhile to extend the method of quantifying and predicting order to non- linear 
models. Non-linear stochastic modeling can be very cumbersome and difficult to implement, but 
some preliminary results are presented here. The same system as in Sec. 3.3.2 and ref. [25] is 
modeled here, and the same parameters are used. 



4.1 ID IMultiscale-lVIultitime Expansion 

First, the results of a ID non-linear model with stochastic evolution is presented. As a first attempt 
at non-linear modeling, two approximations are made. First, the elastic and surface energy parts are 
completely linearized. Second, the wetting potential, W{T-C) is treated using a multiscale-multitime 
expansion. [29, 50] Full details of the model are omitted out of space considerations and because 
these are preliminary results. 

Based on ref. [27], the wetting potential is chosen to be W{H) = 2.314 x lO'^/H erg/cm^ 
with Hin cm. This gives a critical film height of 4 monolayers = 1.132 nm. The simulated film has 
an average film height of = 1.358 nm giving /5 = 0.1447. The simulation cell size is 19.68 fim, 
large enough to hold 513 dots of average size Lq = 2n/ko = 38.4 nm. The multiscale-multitime 
expansion uses an expansion in a scale variable e to create a perturbation-like series. Additionally, 
fast oscillating Fourier components of WCH) are discarded so that the natural length scale is the 
average size of a single dot, Lq . To fourth order in e, one obtains a set of two coupled partial 
differential equations. [29, 50] These equations are solved using Stochastic Euler numerical inte- 
gration [41,42] implemented with Mathematica. [51] Computational efficiency is greatly enhanced 
by the multiscale-multitime expansion, but of course, accuracy and fidelity to the original model 
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Figure 4: Time dependence of number of correlated dots Lcov/Lq vs. dimensionless time, t/tc for 



ID stochastic film evolution (Sec. 4.1 1. Both linear model and multiscale expansion results are 
shown. 



is partially sacrificed. Correlation lengths are calculated from the peak widths of the spectrum 
function (A A;), using Lcor = 1/AA;. The number of dots in a row that form a recognizably periodic 
structure is jj^dots = 2Lcor/ Lq. The time evolution of this number is plotted for both the linear 
model and the stochastic simulation (Fig.|4]). The linear model has a correlation length that grows 
indefinitely. The non-linear model has a correlation length that grows to a peak value and then 
shrinks. In this case, the onset of ripening ruins the SAQD order. The onset of ripening in this 
model corresponds to the "blow-up solution" in ref. [29] 



4.2 2D non-linear deterministic model 

A similar result is obtained for a 2D deterministic non-linear model. This model treats the surface 
energy and wetting energy in full non-linear fashion. The non-local elastic part is found to cubic or- 
der in the film height fluctuation in h via a perturbation series. The stochastic initial conditions are 
sampled white noise with an initial atomic scale roughness, corresponding to = 0.0403 nm^. 
[25] The critical height for the 2D-to-3D-growth transition is Tic = 1132 nm. Two initial aver- 
age film heights are used to investigate the trend predicted in Fig. [s} 7i = 1.257Yc = 1-415 nm 
{(3 = 0.1280) and H = l.VHc = 1.245 nm (/3 = 0.1878). The simulation ceU size corresponds 
to 17 dots squared = 289 dots. The time evolution equations are solved using the native numerical 
differential equation solver in Mathematica. [51] The correlation lengths vs. dimensionless time 
are plotted for both cases in Fig. [5] In both cases, the correlation length increases early on while 
fluctuations are small, reaches a peak value and then decreases due to ripening. The peak value of 
the correlation length is greater for the second case with P closer to the optimal value of 1/4. The 
2D non-linear deterministic model further substantiates the theory that order develops during the 
early growth stages and then is diminished during ripening. Furthermore, the trend predicted by the 
linear order model is it least qualitatively corroborated because growth near the critical threshold 
enhances the peak order of SAQDs according to the 2D non-linear model. 
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Figure 5: Dimensionless correlation length vs. dimensionless time for 2D non-linear deterministic 
model with stochastic initial conditions. Results for two different average film height are reported, 
n = 1.25Kand7^ = 1.1^ 

5 Conclusions 

A linear stochastic model of SAQD order has been presented as an extension of a previous linear 
deterministic model of the order of epitaxial self-assembled quantum dots (SAQDs). In addition, 
some preliminary results from non-linear stochastic and non- linear deterministic models have been 
presented to substantiate the significance of, extend and clarify the linear models of SAQD order. 
The presented numerical examples were based on a very simple SAQD model, and there has been 
much advancement in SAQD growth modeling; however, the presented procedure should apply 
equally well to a wide variety of models with various phenomenological assumptions and help to 
augment them and quantitatively extract order predictions. The current stochastic model should be 
augmented in the future to reflect these advances. Some adaptation of the method ought to apply 
to attempts to engineer SAQD order as well, such as substrate patterning or growing multilayers of 
SAQDs. As with the previous deterministic model, two correlation lengths are found, longitudinal 
L|| and transverse L± . The transverse correlation length appears to be the most limiting, and thus 
should be used to estimate order. It is found that if a wetting potential is incorporated that is suffi- 
ciently soft, growth near the 2D-3D transition critical film height enhances SAQD order; however, 
this enhancement would require very precise experimental control to implement. Nevertheless, it 
demonstrates how the presented methods might apply to other attempts to optimize SAQD growth 
and could help engineer those processes. It was also found that the previous deterministic model 
can be made approximately equivalent to the present stochastic model by choosing the appropri- 
ate initial conditions. Preliminary non-linear modeling appears to corroborate these claims, at 
least qualitatively. A quantitative comparison is still needed. The method to extract SAQD order 
should help with phenomenological model development as the correlation lengths and possibly 
other statistical characterization should facilitate quantitative tuning of phenomenological models 
to experiments. The models presented here apply to the nucleationless mode of SAQD formation; 
however, the inclusion of thermal fluctuations in non-linear models should facilitate a conceptual 
and/or mathematical unification of models of SAQD thermal nucleation and the nucleationless 
mode. 

Thanks to L. Fang for a critical reading of this manuscript. 
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A Derivation of Eq. |7| 

The two-point correlation function in reciprocal space is 

Ckk'(t) = (/ik(t)/^k'W*)- 

Note that at time t = 0, Ckk'(O) = {hk{0)hy,/{0)*) = 0. The time evolution of Ckk(i) can be found 
using the stochastic chain rule (Ito's lemma) and then taking the ensemble average. [41,42] 

StCkk'^ = (^k + ^k')ak'(t)... (14) 
...+ ^^(2nVhT)5^{k-k'). (15) 

The thermal fluctuations only contribute if k = k'. Since initially Ckk'(O) = 0' '^^^ expect 
Ckk'(i) to be non-zero only if k = k'. Thus, the two-point correlation function is determined 
completely by the ensemble averaged spectrum function, Ck(t) as in ref. [25], 

Ckk'W=Ck(t)5'(k-kO. (16) 



This results is only strictly true for the linearized equation. From Eq. 14 the time evolution equation 
of the spectrum function is found by inspection to be Eq. |7] 
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